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ABSTRACT 


Aleutian  Islands  signal  and  noise  data  were  analyzed  to  deter¬ 
mine  the  tele  seismic  recording  capability  of  the  Ocean-Bottom  Seismograph 
(OBS).  From  this  analysis  it  was  concluded  that  both  noise  and  signal  levels 
vary  with  OBS  location,  and  the  lowest  noise  sites  are  in  deep  water  far  from 
land.  The  noise  spectrum  is  sharply  peaked  at  1  Hz;  as  the  frequency  increases 
to  2  Hz,  the  levels  are  down  10  to  20  db.  Thus,  ocean-bottom  seismographs 
can  perceive  higher-frequency  events  up  to  a  full  magnitude  better  than  they  can 
perceive  1 -Hz  events  (important  in  detecting  explosive  events  at  teleseismic 
distances).  OBS  perceptibility  is  quite  variable  at  1  Hz  because  of  the  varia¬ 
bility  of  signal  level  with  location  and  noise  level  with  location  and  time.  It 
appears  that  events  of  magnitude  5.  0  or  greater  can  be  detected  by  at  least  some 
of  the  OBS  stations,  provided  that  the  events  are  40°  or  less  from  the  arrav. 
Events  of  magnitude  6.  0  or  greater  usually  can  be  detected  at  all  epicentral 
distances.  The  noise  field  determines  the  type  of  processing  which  can  be  used. 
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Presented  in  this  report  are  results  of  an  analysis  of  the  Aleu¬ 
tian  Islands  signal  and  noise  data  which  was  conducted  to  determine  the  tele- 
seismic  recording  capability  of  the  Ocean-Bottom  Seismograph  (OBS).  Data 
used  in  the  analysis  were  collected  during  the  3 -phase  1967  Aleutian  Islands 
Ocean- Bottom  Seismographic  Experiment.  A  total  of  349  noise  samples  taken 
from  all  three  phases  were  used  in  the  noise  analysis;  eight  teleseisms  record¬ 
ed  during  Phase  I  were  used  in  the  signal  studies.  Signal  information  con¬ 
tained  in  the  OBS  Preliminary  Analysis  Report  was  also  used. 

Average  noise  levels  at  1  Hz  ranged  from  90  m|i  to  1  (J. ;  the  lower 
noise  levels  were  observed  at  deep  stations  far  from  land.  Except  for  the 
shallow-water  near-shore  sites  where  significant  high-frequency  power  was 
observed,  the  noise  spectrum  peaked  sharply  at  1  Hz;  as  the  frequency  in¬ 
creased  to  2  Hz,  levels  were  down  10  to  20  db  from  the  peak.  Noise  levels 
at  a  site  varied  with  time  due  to  the  effects  of  surface  weather  conditions. 

The  maximum  variation  was  18  db,  with  typical  variations  of  12  db.  Thus, 
the  detection  level  at  a  given  site  varied  as  a  function  of  time  up  to  0.9  mag¬ 
nitude  units.  High-frequency  energy  observed  at  shallow-water  near-shore 
sites  appeared  to  be  related  to  tidal  effects.  Coherence  between  both  the 
pressure  and  vertical  traces  and  the  vertical  and  horizontal  traces  was  low 
over  the  entire  frequency  band,  indicating  that  the  noise  field  was  isotropic 
and  composed  of  several  modes. 

Teleseismic  signal  strength  varied  as  much  as  12  db  (0.  6  mag¬ 
nitude  units)  across  the  200-km  Phase  I  array,  with  6-db  variation  being  typ¬ 
ical.  Some  sites  had  consistently  high  signal  amplitudes;  however,  variations 
appeared  to  depend  partially  on  the  location  of  the  event  as  well. 
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Comparison  of  average  OBS  and  USC&GS  magnitudes  for  five 
events  showed  that  signal  amplitudes  on  the  ocean  bottom  were  slightly  higher 
(0.  2  magnitude  units)  than  on  land.  The  average  period  of  the  tele  seismic 
arrivals  observed  (20°  to  70°  range)  was  0.6  sec  (1.67  Hz),  which  is  some¬ 
what  higher  than  normal.  Because  noise  levels  generally  decrease  rapidly 
above  1  Hz,  higher -frequency  events  are  significantly  better  perceived  by  the 
OBS,  as  probably  reflected  by  the  sample  of  events  used.  Fourteen  of  15  events 
with  a  reported  magnitude  of  5.0  or  better  occurring  between  10°  and  40°  from 
from  the  OBS  stations  were  recorded  by  at  least  one  OBS.  Beyond  40°, 
some  events  were  recorded,  generally  those  which  exhibited  higher  fre¬ 
quency  energy. 


Because  the  noise  field  was  incoherent,  signal  enhancement 
techniques  were  limited  to  summing  and  differencing  the  pressure  and  vertical 
traces  and  to  bandpass  filtering.  About  3-db  signal-to-noise  improvement 
could  be  obtained  using  the  sum-and-difference  techniques,  providing  the  signal 
to-noise  ratios  on  the  pressure  and  vertical  traces  were  similar.  For  the 
higher -frequency  arrivals,  bandpass  filtering  improved  the  signal-to-RMS- 
noise  ratio  and  thus  enhanced  identification  of  the  arrivals. 

The  following  major  conclusions  can  be  drawn. 

•  Both  noise  and  signal  levels  vary  with  OBS  location. 

The  lowest  noise  sites  are  in  deep  water  far  from 
land,  where  noise  levels  of  100  to  200  bat  1  Hz  are 
typical.  Signal  variations  appear  to  be  caused  by  a 
combination  of  site-  and  source -location  effects  and 
are  typically  6  db  across  a  200-km  array. 

•  The  noise  spectrum  is  sharply  peaked  at  1  Hz;  as 
the  frequency  increases  to  2  Hz,  the  levels  are 
down  10  to  20  db.  Thus,  ocean-bottom  seismo¬ 
graphs  can  perceive  higher-frequency  events  up 
to  a  full  magnitude  better  than  they  can  perceive 
1-Hz  events.  Because  explosive  sources  are 
generally  relatively  rich  in  high-frequency  energy, 
the  OBS  would  be  able  to  detect  magnitude  4.  8  to 
5.0  explosive  events  at  tele  seismic  distances. 
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•  Due  to  the  variability  in  both  signal  and  noise 
levels  with  location  (and  with  time  for  the  noise), 
the  OBS  perceptibility  is  quite  variable  at  1  Hz. 

On  the  average,  it  appears  that  events  of  magni¬ 
tude  5.  0  or  greater  can  be  detected  by  at  least 
some  of  the  OBS-array  stations,  provided  that 
the  events  are  40°  or  less  from  the  array. 

Events  of  magnitude  6.0  or  greater  usually  can 
be  detected  at  all  epicentral  distances. 

•  The  type  of  processing  which  can  be  used  depends 
on  the  noise  field.  In  the  worst  case  (incoherent 
noise),  a  signal-to-noise  improvement  of  3  db  can 
often  be  achieved. 
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SECTION  II 


ALEUTIAN  ISLANDS  NOISE  ANALYSIS 

A  detailed  noise  analysis  of  1967  Aleutian  Islands  OBS  data 
was  conducted  to 


•  Define  the  ambient  noise  field  characteristics 
(frequency,  power,  coherence,  etc.  )  within 
the  Aleutian  OBS  array 

•  Relate  the  noise  field  to  the  environment  (water 
depth,  geographic  location,  weather,  and  tides) 

A  total  of  349  3-min  noise  samples,  representing  19  OBS  sta¬ 
tions  and  approximately  1900  km  of  inline  array,  was  digitized  for  this  study 
(Figures  II- 1  and  II-2).  Of  these,  209  were  taken  daily  at  2300  GCT  [1200 
Bering  Daylight  Savings  Time  (BDST)]  and  represent  simultaneous  recordings 
during  each  phase  of  the  operation  (Table  II- 1,  Figure  II - 3) .  The  remaining 
140  samples  were  simultaneous  recordings  from  stations  S2  and  S10  taken 
at  2-hr  intervals  beginning  at  1225  GCT  on  3  September  and  continuing  through 
0625  GCT  on  9  September  1967. 

A.  DESCRIPTION  OF  THE  DATA 

Figure  II-4  outlines  the  order  of  processing  and  analysis  fol¬ 
lowed  for  this  study.  The  data  ensemble  was  selected  after  reviewing  existing 
film  transcriptions  of  the  field  data.  The  0.  0075-ips  recorded  data  were  then 
digitally  sampled  at  16  samples /sec  (0.0625  sec/sample),  resulting  in  an 
8-Hz  folding  frequency.  Antialias  filters  using  a  6-Hz  corner  frequency  and 
a  35-db/octave  rolloff  rate  were  applied  prior  to  digitization. 

Processing  of  these  data  included 

*•  •  Computation  of  power-density  spectra 

•  Computation  of  root-mean- square  (RMS)  values 

•  Visual  period-amplitude  measurements  from  film 
playbacks  with  subsequent  conversion  to  ground  motion 
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ALEUTIAN  NOISE  DATA  ENSEMBLE 
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DATA  DIGITIZATION 


DATA  DISPLAY 


Figure  II-4.  Data  Processing  and  Analysis  Flow  Chart 
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Power-density  spectra  were  computed  for  the  vertical  data 

traces  (xl  and  xlO)  on  each  of  the  daily  samples.  The  spectra  were  computed 

from  slightly  more  than  3  min  of  data  and  have  a  frequency  resolution  of  about 

0.  05  Hz.  To  estimate  the  playback  and  total  system  noise,  spectra  also  were 

computed  on  open-input  recordings  during  transcription  to  digital  format  and 

2 

on  low-level  or  inoperable  component  traces  from  the  field  data.  Figure 
II- 5  shows  that  these  noise  levels  are  substantially  below  seismic  noise  levels 
(shown  later). 


RMS  computations  (a)  were  obtained  for  both  the  daily  and  2-hr 
samples  by  using 


where  N  is  the  number  of  digitial  samples.  These  values  were  used  to  com¬ 
pare  the  total  noise  power  at  each  station  as  a  function  of  water  depth,  weather, 
location,  etc. 


Visual  period-amplitude  measurements  from  film  seismograms 
of  selected  daily  samples  were  converted  to  millimicrons  of  ground  motion. 
Peak  amplitudes  were  measured  at  1  Hz  because  the  noise  peaked  at  this 
frequency  at  nearly  all  stations. 

B.  SUPPLEMENTARY  INFORMATION 

Daily  weather  maps  of  the  North  Pacific  area  were  obtained 
through  the  Naval  Weather  Service  at  Adak,  Alaska.  Barometric  pressures 
at  Amchitka  Island  (recorded  at  0100  local  time)  were  derived  from  these 
charts  and  used  as  weather  indicators  for  the  OBS  arrays. 
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Figure  II-5.  Power-Density  Spectra  of  Total  System  Noise 


The  weather  maps  also  provided  Adak  tidal  information  (time 

O 

and  level  for  high  and  low  tide)  used  in  this  study  because  previous  work 
attributes  significant  microseismic  activity  to  tidal  phenomenon.  Tidal 
characteristics  at  Adak  (Figure  II - 6 )  were  assumed  representative  for  the 
f  mid- Aleutian  area. 

The  weather  and  tide  data  were  combined  with  existing  OBS 
■*  station  information  (water  depth  and  proximity  to  land)  to  determine  the  fac¬ 

tors  which  influence  the  ocean-bottom  noise  levels. 

C.  DATA  ANALYSIS 

The  19  OBS  stations  used  vary  in  water  depth  from  32  to  3840 
fm  (0.  06  to  7.  02  km)  and  in  distance  from  land  from  6.  0  to  1401  km.  Water 
depth  and  distance  from  land  were  the  main  parameters  considered  in  cor¬ 
relating  absolute  noise  amplitudes  with  the  ocean- bottom  environment. 

Noise  spectra  at  most  sites  were  sharply  peaked  at  about  1  Hz 

(Figures  II  —  7  and  II - 8 ) .  Although  these  spectra  have  been  normalized  to 
2 

l(mp/sec)  /Hz  at  1  Hz,  they  have  not  been  corrected  for  system  response, 
so  they  represent  the  frequency  content  seen  on  the  output  records.  At  2  Hz, 
the  noise  levels  are  10  to  20  db  lower. 

Considerable  high-frequency  energy  existed  at  five  sites  (Sll, 
S19,  S25,  S26,  and  S27)  on  half  of  the  samples  (Figure  II-9).  Figure  II- 10 
shows  the  unprocessed  data  at  Sll  during  periods  of  "normal"  and  high- 
frequency  background  noise.  These  five  stations  were  located  in  less  than 
i  750  fm  of  water  and,  except  for  S19,  were  less  than  30  km  from  land  (Table 

II- 1).  Noise  level  at  S19  was  lower  than  at  the  other  sites,  but  a  distinct 
high-frequency  peak  was  still  evident.  Station  S19  was  located  on  top  of  a 
submarine  mountain  where  the  line  of  OBS  stations  crossed  the  Emperor 
Seamounts.  Gains  at  station  S27  (Figure  II - 9)  were  anomalously  low,  which 
accounts  for  the  apparently  low  noise  levels  at  that  shallow-water  near-shore 
site. 
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Figure  II-7.  Power-Density  Spectra  for  Phase  I,  7  September  1967 
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An  attempt  was  made  to  explain  the  high-frequency  energy  on 

3 

the  basis  of  tidal  phenomena,  as  suggested  by  Defant.  High  and  low  tide 
values,  listed  on  daily  weather  charts,  were  plotted  against  time  for  the 
period  spanned  by  this  experiment  (Figure  II-6).  The  RMS  plots  for  shallow 
Phase  III  stations  (shown  later  in  Figure  11-12)  rise  significantly  during  the 
later  portion  of  the  recording  period.  (Tick  marks  on  the  curves  of  Figure 
II- 6  correspond  to  the  time  at  which  the  data  samples  shown  in  Figures  11-12 
and  11-13  were  taken.  )  This  rise  corresponds  with  a  spring-tide  cycle  (period 
of  maximum  tidal  range)  which  begins  about  15-16  July.  In  Phase  II,  large 
increases  in  the  high-frequency  energy  component  (reflected  by  the  sharp 
increases  in  RMS  levels)  at  Sll  (Figure  11-12)  correlate  with  the  early  por¬ 
tion  of  the  spring-tide  cycles  beginning  29-3  0  July  and  12-13  August.  As  in 
Phase  III,  data  samples  are  representative  of  flood  tide.  The  other  Phase  II 
station  (S19)  recording  high-frequency  energy  was  approximately  1030  km 
southwest  of  Sll,  and  the  validity  of  Adak  Island  tide  data  for  such  a  remote 
location  is  questionable.  However,  the  trend  of  the  S19  RMS  plot  is  in  general 
agreement  with  the  respective  portion  of  the  tide  curve  (high  RMS  value  cor¬ 
responding  to  the  latter  portion  of  the  spring-tide  interval  and  a  relative  low 
amplitude  during  neap  tide).  Thus,  it  appears  that  the  high-frequency  energy 
which  is  present  intermittently  at  shallow-water  sites  can  be  related  to  tidal 
activity. 

In  general,  noise-level  variations  correlated  well  with  weather 
for  all  three  phases  of  the  experiment  and  agree  with  past  observations.  Fig¬ 
ure  II- 11  plots  daily  peak  power  levels  in  the  0.  75-  to  1.  25-Hz  band  for  each 
station,  along  with  daily  barometric  pressure  estimates  for  the  area.  The 
effect  of  poor  weather  on  noise  levels  is  quite  clear  for  Phases  I  and  III,  where 
the  pressure  lows  and  noise-level  highs  correlate  very  well.  The  Phase  II 
north  sites  (Sll,  S14,  S15,  and  S16)  also  correlate  with  the  weather;  however, 
the  Phase  II  south  sites  do  not,  probably  because  they  are  very  far  from  Am- 
chitka  where  the  pressure  readings  were  taken. 
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Typical  noise-level  variations  over  the  time  span  of  several 
days  are  12  db  (0.  6  magnitude  units),  with  up  to  18  db  (0.  9  magnitude  units) 
observed  at  some  sites. 

Figures  11-12  and  11-13  plot  the  RMS  levels  as  a  function  of 
time  at  each  site  for  the  vertical  and  pressure  traces,  respectively.  At  the 
deeper  Phase  I  and  Phase  II  sites,  the  plots  are  similar  to  those  of  Figure 
II- 11  and  have  a  similar  shape  for  the  several  instruments  of  the  phases. 
Phase  III  average  power  plots  reflect  the  high-frequency  energy  present  at 
some  sites  and  are  more  variable. 

Phase  I  stations  illustrate  the  relationship  of  noise  levels  to 
water  depth  and  proximity  to  land.  Water  depth  increases  from  Si  to  S7 
then  decreases  slightly  from  S7  to  S10.  The  stations,  of  course,  become 
progressively  farther  from  land.  Noise  levels  generally  decrease  with  in¬ 
creasing  water  depth  and  distance  from  land.  Vertical  data  from  station  S4 
appear  to  be  anomalously  low;  however,  this  probably  is  due  to  instrumental 
gain  problems  because  energy  on  the  pressure  component  has  the  expected 
level.  Pressure  data  from  station  S2  also  appear  low,  but  pressure  data 
from  this  unit  historically  have  been  low,  probably  due  to  its  somewhat 
low  sensitivity.  Vertical  data  from  station  S9  are  more  variable  than  data 
from  other  stations.  Both  the  pressure  and  vertical  components  of  S10  ap¬ 
parently  were  dead  (amplifier  malfunction)  until  6  September,  when  they 
appear  to  have  become  operative,  although  some  problem  may  still  exist 
because  the  vertical  and  pressure  plots  do  not  have  the  same  position  relative 
to  the  plots  from  other  sites. 

Visual  amplitude-period  measurements  from  selected  data 
samples  representing  maximum  and  minimum  1-Hz  energy  levels  at  each 
station  were  converted  to  millimicrons  of  ground  motion  and  are  displayed 
relative  to  water  depth  and  distance  from  land  in  Figure  11-14. 
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Stations  exhibiting  an  undefined  maximum  limit  had  a  large  high- 
frequency  energy  component  during  most  of  the  recording  period  (S26,  S25, 

51  1)  which  obscured  the  low-frequency  energy.  Stations  S28  and  S31  —  having  a 
similar  environment  and  being'  relatively  free  of  high-frequency  energy  — ■  had 
large  1-Hz  energy  levels  (0.  65  to  1.  3 0 (J. ) ,  so  it  is  assumed  that  this  also  holds 
true  for  the  stations  displaying  high-frequency  energy.  The  high-frequency 
energy  at  S19  was  low  enough  to  allow  identification  of  the  1-Hz  energy,  with 
the  resulting  measurement  indicating  low  energy  levels  ( 137-m|i  average  noise 
amplitude)  for  the  1-Hz  ambient  noise.  Ground-motion  average  values,  with 
the  exception  of  shallow-water  near-shore  observations,  were  usually  100  to 
200  mH  which  is  comparable  to  those  observed  from  1966  Kurile  Islands  data 
(Table  II-2).4 

A  detailed  examination  of  RMS  noise  variations  with  time  for 
shallow-  and  deep-water  environments  was  provided  by  2-hr  samples  from 

52  and  S10  data,  respectively.  The  RMS  values  of  these  samples  (Figure  11-15) 
agree  with  the  daily  samples  (tick  marks),  but  they  indicate  that  the  daily  sam¬ 
ples  do  not  necessarily  represent  maximum  and  minimum  noise  levels  for  a 
given  day.  Variations  at  S2  were  12  to  15  db  over  the  time  period  (0.  6  to  0.  7 
magnitude  units);  S10  covered  a  shorter  time,  but  similar  variations  would 
seem  likely,  based  on  the  agreement  between  the  curves.  Station  S10  led 
station  S2  by  2  to  3  hr,  which  suggests  that  the  weather  approached  the  array 
from  the  southwest  at  approximately  25  to  50  knots  during  this  time  period. 

This  observation  is  consistent  with  the  characteristic  weather  patterns  in  the 
area,  as  obtained  from  weather  maps. 

D.  COHERENCE  STUDY 

Stations  S4,  S5,  SI 6,  and  S21  were  used  to  check  the  noise 
coherence  between  components.  Vertical/pressure  coherences  were  com¬ 
puted  for  all  sites;  vertical/horizontal  coherences  were  computed  for  SI 6 
and  S21.  Vertical/pressure  noise  coherences  were  very  low  for  the  Aleutian 
data  over  the  entire  frequency  band  (Figure  11-16).  This  result  contrasts 
sharply  with  the  high  vertical/pressure  coherences  observed  in  the  Gulf  data 
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but  agrees  with  the  results  of  Schneider  et  al  for  Pacific  ocean-bottom  seis¬ 
mic  data  collected  during  1963  in  the  Aleutian  Islands  area.  Vertical/hori¬ 
zontal  coherenc  s  were  also  low,  which  is  characteristic  of  ocean-bottom  noise 
data.  Thus,  it  appears  that  the  noise  field  in  the  Aleutian  area  is  not  only 
isotropic  but  also  is  composed  of  several  modes. 


Table  II-2 

ALEUTIAN  ISLANDS  -  KURILE  ISLANDS 
NOISE-AMPLITUDE  COMPARISON 


Aleutian  Islands 

Kurile  Islands 

Station 

Depth 

(fm) 

Average 

Amp 

(mp.) 

Station 

Depth 

(fm) 

- L - 

Average 

Amp 

(mp) 

26 

32 

366 

- 1 — 

■V 

27 

36 

65 

28 

95 

1045 

4 

230 

136 

25 

101 

780 

4A 

230 

102 

11 

320 

229 

7A 

330 

273 

19 

720 

137 

1 

400 

174 

2 

785 

470 

3A 

600 

102 

31 

862 

352 

1A 

900 

162 

4 

1784 

163 

13 

1320 

92 

16 

2000 

117 

7 

1500 

129 

15 

2025 

169 

11 

1650 

238 

14 

2030 

156 

12 

2850 

95 

5 

2580 

402 

10 

3000 

152 

10 

2690 

169 

5A 

3100 

288 

21 

2840 

91 

9 

2850 

188 

9 

3200 

68 

6 

2980 

299 

18 

3020 

156 

7 

3840 

188 

‘'Instrument  malfunction 
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SECTION  III 


EARTHQUAKE  SIGNAL  ANALYSIS 

A.  DESCRIPTION  OF  THE  DATA 

To  investigate  the  utility  of  certain  signal -enhancement  tech¬ 
niques,  representative  samples  of  teleseismic  P-wave  arrivals  at  several 
Phase  I  instruments  were  digitized  in  TIAC  format  from  the  OBS  field  tapes. 
These  tapes  were  then  converted  to  a  format  compatible  with  the  IBM  System 
360  for  process'. lg.  The  instruments  eelected  were  S2,  S4,  S5,  S6,  S7,  S8, 
S9,  and  S10.  Table  III—  1  shows  the  eight  earthquakes  studied  (designated 
earthquake  1  through  earthquake  8),  their  magnitudes  and  locations,  and  the 
delta  and  azimuth  from  each  earthquake  to  each  station.  The  earthquakes 
span  a  time  period  of  about  3  days. 

A  visual  analysis  of  the  earthquake  recordings  shows  that  sig¬ 
nal  characteristics  vary  across  the  array.  Large  variations  in  signal  ampli¬ 
tude  (for  a  given  earthquake)  are  noted  between  stations  separated  by  as  little 
as  20  km  (Table  HI-2).  Station  S5  tended  to  have  the  largest  amplitudes  and 
S7  the  smallest,  but  the  pattern  varied  from  event  to  event.  This  can  be  seen 
by  comparing  signal  amplitudes  across  the  array  for  earthquakes  3  and  7 
(Figures  III-  1  and  III  -  2) .  The  P-wave  angle  of  incidence  and  the  azimuth  for 
these  two  events  are  quite  different,  suggesting  that  lateral  changes  in  crustal 
and  upper-mantle  structure  in  the  area  also  affect  the  observed  variations. 

Table  III  -  3  compares  the  pressure/vertical  (P/V)  amplitude 
ratios  of  the  noise  and  the  amplitude  ratios  of  the  direct-to- reflected  (water - 
bounce)  signals  on  the  vertical  seismometer.  Theory  predicts  that  an  uncon- 
solidated  bottom  will  give  a  noise  P/V  ratio  of  less  than  unity.  For  this  type 
of  bottom,  the  acoustic  impedance  contrast  between  the  water  layer  and  the 
ocean  bottom  would  be  small,  and  essentially  all  energy  would  be  transmitted 
through  the  interface.  Under  these  conditions,  the  reflected  arrival  is  larger 
than  the  initial  arrival  on  the  vertical  seismometer. 
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Unit  not  down 


Table  III -2 

SIGNAL -AMPLITUDE  VARIATION 


Relative  Amplitude  for  Earthquake 

Station 

2 

3 

5 

6 

7 

S2 

39 

A 

19 

16 

22 

S4 

32 

13 

16 

15 

27 

S5 

34 

26 

28 

30 

33 

S6 

29 

,  19 

B 

B 

33 

S7 

19 

11 

10 

15 

9 

S8 

28 

22 

23 

19 

15 

S9 

25 

A 

14 

C 

C 

S10 

28 

22 

A 

C 

C 

Maximum 

Variation 

2.  1 

2.4 

2.  8 

2.  0 

3.  7 

A  —  Amplitude  too  small  to  measure 

B  —  Event  not  analyzed  because  the  tape  stopped 

C  —  Event  not  analyzed  because  the  seismometer  was  not  on 
on  the  bottom 
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Figure  III—  1 .  Beamsteer  Technique  Applied  to  Earthquake  3,  Averaged 
Time  Differences  Employed  (scale:  1000  TIAC  units/in.) 
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Figure  HI-2.  Beamsteer  Technique  Applied  to  Earthquake  7,  Averaged 
Time  Differences  Employed  (scale:  2000  TIAC  units /in. ) 


Table  III -3 


COMPARISON  OF  P/V  NOISE  RATIOS  AND  DIRECT -TO-REFLECTED 
SIGNAL -AMPLITUDE  RATIOS  ON  THE  VERTICAL  SEISMOMETER 


Station 

P/V  Noise  Ratio 

D/R  Signal- 
Amplitude  Ratio 

S2 

<  1 

<  1 

S4 

>  1 

>  1 

S5 

~  1 

~  1 

S6 

~  1 

~  1 

S7 

>  1 

>  1 

S8 

>1 

>  1 

S9 

<  1 

<1 

S10 

<  1 

<  1 

In  all  cases  where  P/V  <  1,  the  direct-to-reflected  (D/R)  ampli¬ 
tude  ratio  is  also  less  than  unity.  Conversely,  when  P/V  >  1,  D/R  is  also 
>  1.  Figure  HI-2  best  illustrates  this:  S2  has  higher  amplitude  vertical  noise 
and  a  significantly  larger  amplitude  reflected  arrival;  S8  shows  the  opposite 
effect.  From  Table  III-3,  it  therefore  appears  that  the  bottom  was  unconsol¬ 
idated  at  S2,  S 9,  and  S10  and  relatively  consolidated  at  S4,  S7,  and  S8.  Sta¬ 
tions  S5  and  S6  appear  to  fall  between  these  two  sets.  Note  also  that,  for  all 
events,  the  stations  with  a  low  P/V  noise  ratio  tended  to  have  the  best  signal  - 
to-noise  ratio  (especially  S8). 
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Table  III  -  4  compares  USC&GS  magnitudes  with  the  average  mag¬ 
nitudes  computed  using  the  OBS  data.  All  events  are  in  the  20°  to  70°  distance 
range  and  were  recorded  at  more  than  a  single  station.  The  OBS  magnitudes 
averaged  0.  2  magnitude  units  higher.  Data  from  events  less  than  20°  were  ex¬ 
cluded  because  they  showed  considerably  more  scatter  (in  both  directions), 
probably  because  the  phase  being  observed  (i.  e.  ,  travelpath  of  the  first  arrival) 
varied  from  event  to  event.  This  small  sample  suggests  that  teleseismic  signal 
amplitudes  on  the  ocean  bottom  are  about  equal  to  or  slightly  higher  than  thoo e 
at  land  stations. 


Table  III -4 


COMPARISON  OF  OBS  AND  USC&GS  MAGNITUDES 


Event 

Date 

(1967) 

Time 
(h,m,  s) 

OBS 

Station 

Magnitude 

OBS 

Average 

OBS 

USC&GS 

Fiji  Islands 

Aug  1  2 

09:39:44.  3 

S21 

5.  9 

6.0 

5.  8 

SI  9 

5.7 

SI  8 

5.  7 

Sll 

6.  5 

SI  2 

6.  2 

S14 

6.  6 

SI  5 

5.7 

Honshu  Island 

Aug  !  3 

20:06:50.  6 

S21 

5.  9 

6.4 

6.0 

SI  9 

6.  5 

SI  8 

6.  6 

S14 

6.  1 

SI  5 

6.  7 

Kurile  Islands 

Sept  1 

22:42:01.  8 

S5 

5.4 

5.  5 

5.4 

S2 

5.6 

Kermadec 

Sept  7 

11:08:1  3.  2 

S10 

5.  5 

5.  3 

5. 1 

Islands 

S8 

5.  3 

S7 

5.  2 

Mariana 

Sept  9 

08:37:50.4 

S8 

5.  3 

5.  3 

5.  2 

Islands 

S7 

5.2 

S4 

5.1 

S5 

5.  5 

_ 
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B.  DETECTION  RESULTS 


Table  III- 5  is  a  histogram  which  shows  the  number  of  reported 
events  detected  by  the  OBS  arrays  during  each  phase  as  a  function  of  epicentral 
distance.  Between  0°  and  10°,  all  reported  events  were  included;  beyond  10°, 
those  events  with  magnitude  5.  0  and  larger  were  included.  (All  but  four  events 
were  between  magnitudes  5.  0  and  6.  0).  Detections  were  obtained  from  the 
Preliminary  Analysis  Report.  ^  Also,  since  there  may  be  some  reported  events 
which  were  not  positively  associated,  actual  detectability  may  be  slightly  better 
than  indicated. 

Between  0°  and  10°,  all  but  one  reported  event  was  detected; 
between  10°  and  40°,  14  of  15  events  were  detected*  (These  events  generally 
have  shorter  periods  than  events  from  longer  ranges.)  Beyond  40°,  a  few 
scattered  events  (generally  those  deep  events  with  shorter  periods)  were 
detected.  Thus,  it  appears  that,  out  to  about  40°,  events  of  magnitude  5.  0 
and  greater  would  be  detected  by  at  least  some  of  an  OBS  array;  beyond  40°, 
only  very  large  events  or  events  with  relatively  short  periods  would  be  detected. 

C.  SUM  AND  DIFFERENCE  AND  TIME-SHIFTED  SUM  TECHNIQUE 

An  earlier  investigation  of  pressure/vertical  velocity  relation¬ 
ships  for  the  OBS  established  that  formation  of  sum  and  difference  traces  from 

the  pressure  and  vertical  velocity  seismograms  should  permit  the  separation 

6 

of  direct-arrival  energy  from  that  reflected  from  the  water  surface.  The 
pressure  and  vertical  velocity  traces  are  out  of  phase  for  the  direct  P-wave 
arrival  and  in  phase  after  reflection  from  the  water  surface.  Therefore, 
subtracting  the  pressure  trace  from  the  vertical  trace  should  enhance  the 
direct  arrival  and  suppress  the  reflected  arrival;  adding  the  pressure  and 
vertical  traces  should  suppress  the  direct  arrival  and  enhance  the  reflected 
arrival.  Also,  the  signal-to-noise  ratio  for  each  resultant  trace  should  be 
improved  by  3  db  (i.  e.  ,  \fL)  if  the  noise  is  incoherent  between  the  two  input 
channels. 
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Epicentral  distance  measured  from  center  of  array 
All  seismic  activity  was  <5.0  magnitude 


The  sum  and  difference  traces  should  be  nearly  identical;  the 
sum  trace,  however,  will  be  delayed  with  respect  to  the  difference  trace  by 
an  amount  equal  to  the  2~way  traveltime  in  the  water  layer.  If  the  sum  trace 
is  shifted  forward  by  an  amount  equal  to  this  time  delay  and  added  to  the  dif¬ 
ference  trace,  the  result  should  be  approximately  a  fourfold  increase  in  the 
amplitude  of  the  resultant  signal  over  the  original  incoming  signal. 

A  program  was  written  for  the  IBM  System  360  computer  to 
form  the  difference,  sum,  and  time- shifted  sum  traces  for  input  pairs  of 
pressure/vertical  traces  and  to  display  the  result  by  means  of  a  computer¬ 
generated  plot.  The  amplitude  of  the  time-shifted  and  summed  trace  is 
divided  by  four  before  plotting  so  that  it  may  be  compared  more  readily  with 
the  input  traces. 

Figures  III- 3  and  III  -  4  show  two  such  plots  for  earthquake  7, 
recorded  at  stations  S5  and  S6.  There  is  a  significant  improvement  in  the 
signal-to-noise  ratio  for  each  of  the  three  output  traces  over  the  two  input 
traces.  However,  summing  the  input  traces  does  not  completely  cancel  the 
portion  of  the  incoming  signal  ahead  of  the  water- surface  reflection.  When 
the  sum  trace  is  shifted  and  added  to  the  difference  trace,  the  incompletely 
cancelled  portion  shows  up  as  a  precursor  to  the  incoming  signal,  making  the 
exact  point  of  signal  onset  less  distinct  than  it  is  on  the  difference  trace. 
Similar  comments  apply  to  Figures  III- 5  and  III- 6  which  show  arrivals  from 
earthquake  5  at  stations  S5  and  S9. 

Figures  III-7  and  III -8  show  the  technique  applied  to  signal 
arrivals  from  earthquake  3  at  stations  S7  and  S10.  In  this  instance,  little 
improvement  in  signal-to-noise  ratio  of  the  output  traces  over  the  better  of 
the  two  input  traces  is  achieved  because,  in  each  case,  one  of  the  input  traces 
has  a  significantly  higher  noise  level  than  the  other  and  a  relatively  low  signal 
level. 
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Figure  III-3.  Sum  and  Difference  Technique  Applied  to 
Arrivals  from  Earthquake  7  at  Station  35 
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Figure  III  -  7 .  Sum  and  Difference  Technique  Applied  to 
Arrivals  from  Earthquake  3  at  Station  S7 


Figure  III  -  8 .  Sum  and  Difference  Technique  Applied  to 
Arrivals  from  Earthquake  3  at  Station  S10 
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The  method  was  unsuccessfully  applied  to  earthquake  1,  a 
relatively  weak  event  which  was  not  detected  on  any  of  the  selected  instru¬ 
ments  . 


Figure  III  -  9  shows  the  method  applied  to  the  signal  arrival 
from  earthquake  6  at  station  S7.  In  this  case,  the  input  data  better  satisfy 
the  theoretical  assumptions  of  the  method;  the  input  traces  have  approximately 
the  same  amplitude  and  signal-to-noise  ratio,  and  the  initial  arrival  is  almost 
completely  cancelled  in  the  difference  trace.  The  signal  arrival  on  the  final 
output  trace  shows  a  sharp  onset,  and  the  signal-to-noise  ratio  for  this  trace 
is  between  3  and  6  db  higher  than  it  is  on  either  input  trace. 

In  summary,  this  technique  will  give  at  least  3-db  signal-to- 
noise  improvement  if  the  signal-to-noise  ratio  on  the  pressure  and  vertical 
components  are  similar  and  the  noise  is  uncorrelated  between  the  two  channels. 
In  most  cases,  the  difference  trace  is  better  than  the  other  combinations  be¬ 
cause  no  precursor  is  introduced  to  the  output  trace  (which  is  caused  by  im¬ 
perfect  cancellation  of  the  energy  between  the  initial  arrival  and  water  bounce). 
Occasionally,  however,  the  sum  of  the  difference  trace  and  tirr  e- shifted  sum 
trace  can  be  used. 

D.  BEAMSTEERING 

The  beamsteer  technique  is  a  method  of  aligning  and  summing 
the  seismograms  from  an  array  of  seismographs  to  increase  event-detection 
capability.  If  the  noise  is  incoherent  between  the  input  traces  and  the  wave¬ 
form  is  preserved  across  the  array,  the  method  should  result  in  i/n  improve¬ 
ment  of  the  signal-to-noise  ratio  of  the  output  trace  over  the  input  traces, 
where  n  ir.  the  number  of  input  traces  summed.  Because  the  data  were  readily 
available,  the  method  was  applied  to  the  earthquake  recordings  even  though 
it  is  known  that  array  dimensions  this  large  (160  km)  can  be  expected  to  in¬ 
troduce  critical  signal-alignment  problems,  as  well  as  degradation  due  to 
changes  in  signal  character. 
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Figure  III- 9.  Sum  and  Difference  Technique  Applied  to 
Arrivals  from  Earthquake  6  at  Station  ST 
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The  theoretical  traveltimes  from  a  strong  earthquake  to  each 
array  station  were  computed;  from  these,  the  theoretical  arrival  times  of  the 
signal  at  each  unit  were  determined.  The  difference  between  the  theoretical 
and  true  arrival  time  was  noted  for  each  station.  Time  differences  are  as¬ 
sumed  to  be  due  to  crustal  variations  or  other  localized  effects  which  are 
characteristic  of  the  station  and  independent  of  the  earthquake  parameters. 

In  this  investigation,  the  time  differences  at  each  station  for  all 
the  visually  detectable  earthquakes  —  earthquakes  2,  3,  5,  6,  7,  and  8  —  were 
averaged  to  determine  the  time-difference  ( static- correction)  characteristic 
of  each  station.  Static- correction  information  is  input  to  a  computer  program 
which  shifts  the  traces,  sums  them,  and  plots  the  aligned  input  traces  and  the 
sum  trace.  The  pressure  and  vertical  velocity  channels  are  used  as  dual  in¬ 
puts  for  each  station,  the  pressure  traces  being  inverted  during  summation  to 
account  for  the  180°  phase  difference  between  :he  pressure  and  vertical  traces. 
Also,  to  facilitate  comparison  with  the  input  traces,  ■uhe  amplitude  of  the  output 
trace  is  divided  by  the  number  of  input  traces  before  it  is  plotted. 

Figures  III- 1  and  III  -  2  and  III- 10  through  III- 15  show  the  individ¬ 
ual  seismometers  and  the  beamsteered  outputs  of  this  program  for  earthquakes 
1,  3,  4,  5,  6,  7,  and  8  using  the  static  correction  just  described.  The  method 
did  not  result  in  significant  signal-to-nois e- ratio  improvement  or  onset-detec¬ 
tion  capability  over  the  better  input  traces  due  to  the  following  factors. 

There  is  a  substantial  variation  in  signal  strength  and  signal- 
to-noise  ratio  over  the  array,  and  the  sum  trace  is  degraded  by  the  traces 
with  poor  signal-to-noise  ratios.  A  technique  of  weighting  the  input  traces 
based  on  their  observed  signal-to-noise-ratios  (diversity  stacking)  might  have 
improved  results.  For  example,  Figures  III- 12  and  III- 13  show  the  program 
output  for  earthquake  5  with  and  without  the  vertical  trace  of  unit  S8 .  In  this 
case,  the  output  trace  is  virtually  unaltered  when  the  S8  vertical  trace  is 
omitted,  even  though  this  trace  has  a  significantly  better  signal-to-noise  ratio 
than  any  of  the  others. 
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Figure  III- 11.  Beamsteer  Technique  Applied  to  Earthquake  4,  Averaged 

Time  Differences  Employed  (scale:  1000  TIAC  units /in. ) 
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Figure  III- 12.  Beamsteer  Technique  Applied  to  Earthquake  5,  Averaged 

Time  Differences  Employed  (scale:  1000  TIAC  units /in.  ) 
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Figure  III- 13.  Beamsteer  Technique  Applied  to  Earthquake  5,  Averaged 

Time  Differences  Employed,  Vertical  Trace  of  Station 
S8  Omitted  (scale:  1000  TIAC  units/in.  ) 


III- 21 


science  services  division 


"d  . 
«  a 

OU  .!-( 
nJ 

^  m 
<U  +J 

>  3 

C  3 


III-22 


science  services  division 


%T 

• ,  t 

-  -TffnwrtiwtfMnon  mi . „ 

■■■MUMP 

:-SS®f.5K!«!SKasSi»ii 

MernKS1®?, 

Bl 

^«wtiR?BssSi3?8s?sa 

vt 

ii  ‘  ’fe.7W,>i.,FRV51  •  J'  k  I 

Beamsteer  Technique  Applied  to  Earthquake  8,  Averaged 
Time  Differences  Employed  (scale:  1000  TIAC  units /in. ) 


The  sampling  rate  of  the  digitized  data  makes  the  amount  of 

^  . 

the  time  shifts  critical.  For  the  sample  rate  used,  a  difference  of  only  eight 
data  blocks  is  equivalent  to  a  phase  difference  of  180°  at  1  Hz. 

The  average  static  corrections  were  not  sufficient.  Figures 
III- 16  and  III- 17  show  the  plotted  program  outputs  for  earthquakes  3  and  7, 
respectively;  the  traces  are  aligned  and  summed  using  the  time  shifts  indi¬ 
cated  by  the  individual  earthquake  rather  than  the  averaged  characteristic 
time  difference.  In  each  case,  the  sum  trace  using  trace-alignment  informa¬ 
tion  from  the  individual  earthquake  shows  a  higher  amplitude  and  signal-to- 
noise- ratio  than  the  sum  trace  computed  using  the  information  averaged  from 
several  earthquakes.  Thus,  it  appears  that  the  static  corrections  are  a  func¬ 
tion  of  earthquake  location  also  (not  a  surprising  result). 

E.  ZERO-PHASE  BANDPASS  FILTERING 

Presented  here  are  results  of  calculations  aimed  at  improving 
signal -to -noise  ratios  by  simple  bandpass  filtering.  Instrument  response  can 
be  changed  to  afford  some  noise  rejection  (thereby  lowering  the  signal-detection 
threshold  of  the  OBS),  provided  the  signal  and  noise  spectra  peak  at  different 
frequencies.  Change  of  the  instrument  response  was  simulated  in  the  computer 
by  applying  zero-phase  bandpass  filters  to  the  pressure  and  vertical  velocity 
recordings  of  three  earthquake  events. 

Power  spectra  were  computed  for  signals  on  both  the  pressure 
and  vertical  channels  of  S4,  S5,  and  S6  for  earthquake  7  (magnitude  5.8).  That 
event  was  chosen  because  its  high  signal-to-noi  se  ratio  allowed  a  good  esti¬ 
mate  of  the  signal  spectra.  In  addition,  noise  power  spectra  were  computed 
for  noise  samples  immediately  preceding  each  signal.  Calculations  were 
performed  on  the  IBM  System  360  computer  using  the  Cooley- Tukey  algorithm 
(fast  Fourier  transform),  and  CalComp  plots  of  the  power  spectra  were  obtained. 
The  number  of  points  in  each  time  trace  was  1024  and  the  sample  rate  (At)  was 
0.  0625  sec,  giving  a  frequency  increment  (Af)  of  0.  015625  Hz  and  a  gate  length 
of  64  sec. 
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Figure  III- 16.  Beamsteer  Technique  Applied  to  Earthquake  3,  Optimum 

Visual  Shift  Employed  (scale:  1000  TIAC  units/in.) 
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Figure  III- 17.  Beamsteer  Technique  Applied  to  Earthquake  7,  Optimum 

Visual  Shift  Employed  (scale:  2000  TIAC  units/in.) 


The  corresponding  signal  and  noise  spectra  are  plotted  together 
in  Figure  III- 18  to  allow  compar'son  of  signal-to-noise  ratios  as  a  function  of 
frequency.  An  inspection  of  the  plots  shows  that  the  signal  spectra  for  this 
earthquake  peaked  at  about  1.8  Hz.  At  this  frequency,  the  noise  level  on  the 
vertical  component  is  down  between  10  and  15  db  from  the  1-Hz  microseismic 
peak.  If  the  earthquake  had  been  a  full  magnitude  unit  lower  (4.8)  with  the 
same  spectrum,  a  signal-to-noise  ratio  of  10  db  (3:1)  at  1.8  Hz  would  still 
have  existed. 

Table  III  -  6  shows  that  the  average  period  of  teleseismic  arrivals 
was  0.  6  sec  (1.7  Hz).  This  relatively  high  frequency  is  probably  the  resuT;  of 
the  lower  frequency  arrivals  (around  1  Hz)  being  obscured  by  noise  (consider¬ 
ably  higher  at  1  Hz)  so  that  only  the  higher  frequency  events  are  observed.  It 
should  be  pointed  out  that  explosive  sources  generally  have  a  spectrum  which 
is  relatively  rich  in  high-frequency  energy;  therefore,  the  detection  capability 
of  the  OBS  would  be  significantly  better  for  explosions  than  for  earthquakes. 


Table  III -6 

AVERAGE  PERIOD  OF  TELESEISMIC  ARRIVALS 


Delta 

(°) 

Number  of 

Station  Observations 

Average  Period 
( s  ec) 

10-20 

15 

0.21 

20-30 

9 

0.  64 

30-40 

4 

0.  55 

40-50 

4 

0.  55 

50-60 

— 

— 

60-70 

5 

0.  52 

70-80 

6 

0.  67 

T  otal 

28 

0.  60 

(20°  to  80°) 
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Five  zero-phase  shift  digital  filter  operators  (2- sec  in  lenght) 
were  designed.  The  response  of  each  filter  is  described  by  the  cutoff  fre¬ 
quencies  f^,  f  i ,  f^,  and  f^  given  in  Table  III  -  7 ,  where  the  desired  frequency 
response  H(f)  is 


H(f) 

H(f) 

H(f) 

H(f) 

H(f) 


for  0  £  f  £  f 


2 


2 


+  —  COS  TT 


-  -  COS  TT 


f-fl 

for  f 

o 

sf, 

£o‘£l 

o 

1 

for  f 

£f  £f 

f-f3 

1 

for  f_ 

2 

sif  if 

f2  ‘  f3 

2 

3 

r-H 

VI 

VlH 

VI 

for  f^ 

2  At 

Table  III -7 

FILTER  CUTOFF  FREQUENCIES 


Filter 

Frequency  (Hz) 

f 

o 

£i 

f2 

£3 

1 

0.  7 

1.  0 

1.3 

1.7 

2 

0.  9 

1.2 

1.  6 

2.2 

3 

1.  0 

1.6 

2.  0 

2.  6 

4 

1.4 

2.  0 

2.4 

3.  0 

5 

0.  0 

0.  0 

2.  0 

3.  0 
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Filters  3  and  4  were  applied  by  convolution  to  the  digital  re¬ 
cordings  of  earthquake  3.  The  input  traces  are  shown  in  Figure  III- 19  (where 
the  signals  are  vertically  aligned);  the  filtered  outputs  of  filters  3  and  4  are 
shown  in  Figures  III-20  and  III-21,  respectively.  Filter  4  gave  noticeable 
signal- to -noise  -  ratio  improvement  on  some  of  the  noisier  channels  because 
the  filter  passband  was  centered  near  the  signal  frequency  and  the  dominant 
noise  frequency  was  outside  the  passband.  The  ringing  appearance  of  the 
traces  is  due  to  the  narrow  passband  specified  in  the  filter  response. 

Filters  1  through  4  were  applied  to  earthquake  5.  The  input 
traces  are  shown  in  Figure  III-22,  and  the  filtered  outputs  appear  in  Figures 
III  -  2  3  through  III- 26 .  For  this  earthquake  filter  1  was  centered  near  the 
signal  frequency,  which  should  give  the  best  signal-to-noise  improvement. 
However,  the  signal  frequency  was  closer  to  the  dominant  noise  frequency,  so 
the  improvement  was  less  than  that  obtained  for  earthquake  3.  This  filter  did 
reject  a  significant  amount  of  the  low-  and  high-frequency  noise  energy  out¬ 
side  the  passband.  The  low-amplitude  levels  seen  in  Figures  III-25  and  111-26 
for  filters  3  and  4  indicate  low  signal  and  noise  energy  in  those  frequency  bands. 

If  the  OBS  were  located  in  a  highly  seismic  region,  a  large 
number  of  teleseismic  events  would  be  obscured  by  local  and  near-regional 
earthquake  activity.  For  example,  a  total  of  2230  unassociated  station  events 
were  recorded  during  the  1967  Aleutian  Islands  Experiment.  These  local  and 
near- regional  events  would  be  considered  as  undesirable  noise  in  terms  of  the 
capability  to  record  teleseisms.  Generally,  however  (since  their  energy  is 
predominately  high  frequency),  events  could  be  removed  from  the  recordings 
by  changing  the  instrument  response. 

-f 


III- 30 


science  services  division 


Figure  111-27  shows  earthquake  2  masked  by  a  near- regional 
event.  The  predicted  arrival  time  for  the  teleseismic  P-wave  is  indicated 
for  each  station,  the  traces  are  not  aligned.  The  telesismic  P-wave  is  ob¬ 
scured  at  most  stations  by  the  high-frequency  ^ear-regional  event.  Filter  5 
was  applied  to  the  traces,  and  the  output  is  shown  in  Figure  III -28  where  the 
teleseismic  signal  is  aligned  vertically.  The  high-cut  filter  adequately  re¬ 
moved  the  near-regional  event  from  tk-  recordings,  allowing  the  telesiesmic 
P-wave  to  be  observed. 
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Figure  III- 1 9-  Earthquake  3,  No  Filtering  (scale:  1000  TIAC  units/in.) 
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Figure  III-20.  Filtered  Output  of  Earthquake  3,  Filter  3  (scale:  1000  TIAC  units/in.  ) 
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Figure  HI-21.  Filtered  Output  of  Earthquake  3,  Filter  4  (scale:  1000  TIAC  units/in. 
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Figure  III- 22.  Earthquake  5,  No  Filtering  (scale:  1000  TIAC  units/in.) 
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Figure  III-23.  Filtered  Output  of  Earthquake  5,  Filter  1  (scale:  1000  TIAC  units/in. 
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Figure  III-24.  Filtered  Output  of  Earthquake  5,  Filter  2  (scale:  1000  TIAC  units /in. 


1 


III- 39 


science 


I  iiir* 


Figure  HI-26.  Filtered  Output  of  Earthquake  5,  Filter  4  (scale:  1000  TIAC  units/in.) 
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igure  III- 2 3. 


Filtered  Output  of  Earthquake  2,  Filter  5,  Traces 
Vertically  Aligned  (scale:  1000  TIAC  units/in.) 
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APPENDIX  A 


RELATIONSHIP  BETWEEN  OBS  INSTRUMENT 
MAGNIFICATION  AND  TIAC  UNITS 

Recent  ocean-bottom  seismograph  shaketable  tests  were  used 

* 

to  relate  film  amplitudes  to  ground  motion.  The  shaketable  data  were  used 
to  derive  a  similar  relationship  for  the  digitized  data  (TIAC  units/mp). 

Unit  21  shaketable  data  from  0.6  to  12  Hz  were  selected  and 
digitized  using  a  0.  03125-sec  sample  rate.  The  data  were  demultiplexed  and 
transcribed  to  the  IBM  System  360  for  input  to  a  program  which  determined 

•  Mean  (which  was  removed) 

•  Peak  a  mplitudes 

•  Average  peak  amplitude 

•  Average  frequency 

Amplitudes  (Table  A-l)  were  averaged  over  the  20  largest  0- 
to-peak  values  at  a  given  frequency  after  removal  of  the  mean.  Using  the 
known  table  amplitude,  a  conversion  factor  (TIAC  units /  mp)  was  calculated 
and  plotted  (Figure  A-l).  The  xl  and  xlO  channels  are  separated  by  the  ex¬ 
pected  factor  of  about  10  from  1.  5  to  8.  5  Hz.  However,  outside  this  band, 
system  noise  masks  the  calibration  signal  on  the  xl  channel,  so  the  xlO  data 
were  used  to  obtain  the  conversion. 

Using  the  TIAC  units/mp  ratio  calculated  at  1  Hz,  the  power 

2 

spectra  were  converted  to  absolute  units  -  db  relative  to  l(mp)  /Hz  at  1  Hz. 
The  spectra  shown  throughout  the  report  are  in  absolute  units. 


Texas  Instruments  Incorporated,  1968:  Ocean-Bottom  Seismograph  System 
Response  Study,  Aleutian  Islands  Experiment,  Contract  F33657-67-C-1  341 , 
31  Jan. 
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Table  A-l 


UNIT  21  SHAKETABLE  TESTS 


(Vertical  Component  at  -60  db) 


F  requency 
(Hz) 

Table  Magnitude  (p) 

Amplitude  (0-p) 
of  TIAC  Units 

P-P 

°-p 

xl 

xlO 

0.6 

9.6 

4.  8 

48.  038 

78.  270 

0.  7 

9.6 

4.  8 

52. 973 

98.466 

0.  8 

9.6 

4.  8 

44. 039 

119. 486 

0.  9 

9.2 

4.  6 

48.  045 

190. 705 

1.0 

9.0 

4.  5 

61.  704 

274. 469 

1. 2 

9.0 

4.  5 

74. 074 

455. 524 

1.4 

8.  6 

4.  3 

78. 808 

633. 805 

1.  8 

8.4 

4.  2 

97.484 

85  7.  575 

2.  0 

8.4 

4.  2 

106.200 

948.340 

2.  2 

8.  2 

4.  1 

111. 270 

1037.  794 

2.4 

8.  0 

4.  0 

116. 511 

1065. 146 

2.  8 

7.8 

3.  9 

125. 387 

1179. 009 

3.0 

7.  6 

3.  3 

135. 379 

1214. 195 

3.  2 

7.  5 

3.  75 

138.627 

1233. 574 

3.4 

7.  2 

3.6 

122. 720 

1219. 900 

3.8 

7.  0 

3.  5 

133. 847 

1235. 298 

4.  0 

7.  0 

3.  5 

l$0. 324 

1244. 846 

4.  2 

6.  8 

3.4 

132. 954 

1235.414 

4.  4 

6.  7 

3.35 

135.  418 

1202. 243 

4.  8 

6 . 6 

3.3 

135. 204 

1152. 118 

5.  0 

6.4 

3.  2 

130. 300 

1147.  197 

5.  5 

6.  1 

3. 05 

130. 836 

1127.  754 

6.  0 

5.  8 

2.  9 

103. 250 

1015.495 

7.  0 

5.  2 

2.  6 

110. 028 

974. 777 

7.  5 

5.  0 

2.  5 

110. 800 

941. 915 

8.  0 

4.  8 

2.  4 

95. 875 

610. 299 

8.  5 

4.  5 

2.  25 

72. 269 

558. 598 

10.  0 

3.  8 

1.9 

68.  207 

368.  958 

12.  0 

3.2 

1.6 

55. 900 

129.  718 

Amplitude  Normalized 
to  1  mp  and  0  db 
(units  /mp) 


xl 


10. 008 
11. 036 
9.  175 
10.445 
13. 712 
16.461 
18. 327 
23.210 
25. 286 
27.  139 
29. 128 
32. 151 
35.626 
36.967 
34. 089 
38. 242 
37.  235 
39. 104 
40.  42 
40. 971 
40. 719 

42. 897 
35.  603 
42.318 
44. 320 
39. 948 
32. 1 19 

35. 898 
34. 937 


xlO 


16. 306 
20.  514 
24. 893 
41. 458 
60. 993 
101. 228 
147.  397 
204. 185 
225.  795 
253. 120 
266. 287 
302.310 
319. 525 
328.  953 
338.  861 
352. 942 
355.670 
363.357 
358. 879 
349.  127 
358.  499 
369. 755 
350. 1 71 
374. 914 
3  76.  766 
254.  291 
248.  266 
194. 188 
8 1 . 0  74 
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APPENDIX  B 


REPORTS  PREPARED  UNDER  OCEAN-BOTTOM 
SEISMOGRAPHIC  CONTRACTS 


The  following  reports  were  written  under  Contract  F33657-67- 

C - 1341: 


•  Operations  Report,  Aleutian  Islands  Experiment, 
Ocean-Bottom  Seismographic  Experiments,  31 
January  1968 

•  Ocean-Bottom  Seismograph,  System-Response 
Study,  31  January  1968 

•  Preliminary  Analysis  Report,  Aleutian  Islands 
Experiment,  Ocean-Bottom  Seismographic  Ex¬ 
periments,  31  January  1968 

•  Final  Analysis  Report,  Aleutian  Islands  Experi¬ 
ment,  Ocean-Bottom  Seismographic  Experiments, 

31  July  1968 

•  Aleutian  Mantle  Study,  Aleutian  Islands  Experi¬ 
ment,  Ocean-Bottom  Seismographic  Experiments, 

31  July  1968 

•  Noise  Analysis  Report,  Aleutian  Islands  Experi¬ 
ment,  Ocean-Bottom  Seismographic  Experiments, 

31  October  1968 

Reports  written  under  Contract  F33657-68-C- 0242  are  as 

follows: 

•  Final  Report,  Ocean-Bottom  Seismograph  Pro¬ 
duction  and  Gulf  of  Mexico  Data  Analysis,  31 
July  1968 

•  Operation  Manual  for  30-Day  Ocean-Bottom  Seismo¬ 
graph,  August  1968 

Written  under  Contract  F33657-68-C-0875  was  this  report: 

•  Final  Report,  1968  Aleutian  Islands  Experiment, 
Ocean-Bottom  Seismographic  Experiments,  23 
August  1968 
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